NUPARM and NUCGEN: software for analysis and generation of sequence dependent nucleic acid structures

Software packages NUPARM and NUCGEN, are described, which can be used to understand sequence directed structural variations in nucleic acids, by analysis and generation of non-uniform structures. A set of local inter basepair parameters (viz. tilt, roll, twist, shift, slide and rise) have been defined, which use geometry and coordinates of two successive basepairs only and can be used to generate polymeric structures with varying geometries for each of the 16 possible dinucleotide steps. Intra basepair parameters, propeller, buckle, opening and the C6...C8 distance can also be varied, if required, while the sugar phosphate backbone atoms are fixed in some standard conformation ill each of the nucleotides. NUPARM can be used to analyse both DNA and RNA structures, with single as well as double stranded helices. The NUCGEN software generates double helical models with the backbone fixed in B-form DNA, but with appropriate modifications in the input data, it can also generate A-form DNA ar rd RNA duplex structures.

The carboxyl-terminal domain of the p53 protein regulates sequence-specific DNA binding through its nonspecific nucleic acid-binding activity.

The murine p53 protein contains two nucleic acid-binding sites, a sequence-specific DNA-binding region localized between amino acid residues 102-290 and a nucleic acid-binding site without sequence specificity that has been localized to residues 364-390. Alternative splicing of mRNA generates two forms of this p53 protein. The normal, or majority, splice form (NSp53) retains its carboxyl-terminal sequence-nonspecific nucleic acid-binding site, which can negatively regulate the sequence-specific DNA-binding site. The alternative splice form of p53 (ASp53) replaces amino acid residues 364-390 with 17 different amino acids. This protein fails to bind nucleic acids nonspecifically and is constitutive for sequence-specific DNA binding. Thus, the binding of nucleic acids at the carboxyl terminus regulates sequence-specific DNA binding by p53. The implications of these findings for the activation of p53 transcriptional activity following DNA damage are discussed.

Nucleotide sequence of 5S rDNA and localization of the ribosomal RNA genes to metaphase chromosomes of the Tilapiine cichlid fish, Oreochromis niloticus

In this study, we report the cloning and nucleotide sequence of PCR-generated 5S rDNA from the Tilapiine cichlid fish, Oreochromis niloticus. Two types of 5S rDNA were detected that differed by insertions and/or deletions and base substitutions within the non-transcribed spacer (NTS). Two 5S rDNA loci were observed by fluorescent in situ hybridization (FISH) in metaphase spreads of tilapia chromosomes. FISH using an 18S rDNA probe and silver nitrate sequential staining of 5S-FISH slides showed three 18S rDNA loci that are not syntenic to the 5S rDNA loci.

Sp1 up-regulates cAMP-response-element-binding protein expression during retinoic acid-induced mucous differentiation of normal human bronchial epithelial cells.

CREB [CRE (cAMP-response element)-binding protein] is an important transcription factor that is differentially regulated in cells of various types. We recently reported that RA (retinoic acid) rapidly activates CREB without using RARs (RA receptors) or RXRs (retinoid X receptors) in NHTBE cells (normal human tracheobronchial epithelial cells). However, little is known about the role of RA in the physiological regulation of CREB expression in the early mucous differentiation of NHTBE cells. In the present study, we report that RA up-regulates CREB gene expression and that, using 5'-serial deletion promoter analysis and mutagenesis analyses, two Sp1 (specificity protein 1)-binding sites located at nt -217 and -150, which flank the transcription initiation site, are essential for RA induction of CREB gene transcription. Furthermore, we found that CREs located at nt -119 and -98 contributed to basal promoter activity. Interestingly, RA also up-regulated Sp1 in a time- and dose-dependent manner. Knockdown of endogenous Sp1 using siRNA (small interfering RNA) decreased RA-induced CREB gene expression. However, the converse was not true: knockdown of CREB using CREB siRNA did not affect RA-induced Sp1 gene expression. We conclude that RA up-regulates CREB gene expression during the early stage of NHTBE cell differentiation and that RA-inducible Sp1 plays a major role in up-regulating human CREB gene expression. This result implies that co-operation of these two transcription factors plays a crucial role in mediating early events of normal mucous cell differentiation of bronchial epithelial cells.

Nucleotide sequences of an Australian and a Canadian isolate of potato leafroll luteovirus and their relationships with two European isolates

The genomes of an Australian and a Canadian isolate of potato leafroll virus have been cloned and sequenced. The sequences of both isolates are similar (about 93%), but the Canadian isolate (PLRV-C) is more closely related (about 98% identity) to a Scottish (PLRV-S) and a Dutch isolate (PLRV-N) than to the Australian isolate (PLRV-A). The 5'-terminal 18 nucleotide residues of PLRV-C, PLRV-A, PLRV-N and beet western yellows virus have 17 residues in common. In contrast, PLRV-S shows no obvious similarity in this region. PLRV-A and PLRV-C genomic sequences have localized regions of marked diversity, in particular a 600 nucleotide residue sequence in the polymerase gene. These data provide a world-wide perspective on the molecular biology of PLRV strains and their comparison with other luteoviruses and related RNA plant viruses suggests that there are two major subgroups in the plant luteoviruses.

Chiropteran types I and II interferon genes inferred from genome sequencing traces by a statistical gene-family assembler.

BACKGROUND: The rate of emergence of human pathogens is steadily increasing; most of these novel agents originate in wildlife. Bats, remarkably, are the natural reservoirs of many of the most pathogenic viruses in humans. There are two bat genome projects currently underway, a circumstance that promises to speed the discovery host factors important in the coevolution of bats with their viruses. These genomes, however, are not yet assembled and one of them will provide only low coverage, making the inference of most genes of immunological interest error-prone. Many more wildlife genome projects are underway and intend to provide only shallow coverage. RESULTS: We have developed a statistical method for the assembly of gene families from partial genomes. The method takes full advantage of the quality scores generated by base-calling software, incorporating them into a complete probabilistic error model, to overcome the limitation inherent in the inference of gene family members from partial sequence information. We validated the method by inferring the human IFNA genes from the genome trace archives, and used it to infer 61 type-I interferon genes, and single type-II interferon genes in the bats Pteropus vampyrus and Myotis lucifugus. We confirmed our inferences by direct cloning and sequencing of IFNA, IFNB, IFND, and IFNK in P. vampyrus, and by demonstrating transcription of some of the inferred genes by known interferon-inducing stimuli. CONCLUSION: The statistical trace assembler described here provides a reliable method for extracting information from the many available and forthcoming partial or shallow genome sequencing projects, thereby facilitating the study of a wider variety of organisms with ecological and biomedical significance to humans than would otherwise be possible.

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing.

BACKGROUND: There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important. RESULTS: Here we provide a systematic exploration of how well RNA-Seq can identify human coding variants by comparing variants identified through high coverage whole-genome sequencing to those identified by high coverage RNA-Seq in the same individual. This comparison allowed us to directly evaluate the sensitivity and specificity of RNA-Seq in identifying coding variants, and to evaluate how key parameters such as the degree of coverage and the expression levels of genes interact to influence performance. We find that although only 40% of exonic variants identified by whole genome sequencing were captured using RNA-Seq; this number rose to 81% when concentrating on genes known to be well-expressed in the source tissue. We also find that a high false positive rate can be problematic when working with RNA-Seq data, especially at higher levels of coverage. CONCLUSIONS: We conclude that as long as a tissue relevant to the trait under study is available and suitable quality control screens are implemented, RNA-Seq is a fast and inexpensive alternative approach for finding coding variants in genes with sufficiently high expression levels.

Expansion of the alpha 2-adrenergic receptor family: cloning and characterization of a human alpha 2-adrenergic receptor subtype, the gene for which is located on chromosome 2.

Pharmacologic, biochemical, and genetic analyses have demonstrated the existence of multiple alpha 2-adrenergic receptor (alpha 2AR) subtypes. We have cloned a human alpha 2AR by using the polymerase chain reaction with oligonucleotide primers homologous to conserved regions of the previously cloned alpha 2ARs, the genes for which are located on human chromosomes 4 (C4) and 10 (C10). The deduced amino acid sequence encodes a protein of 450 amino acids whose putative topology is similar to that of the family of guanine nucleotide-binding protein-coupled receptors, but whose structure most closely resembles that of the alpha 2ARs. Competition curve analysis of the binding properties of the receptor expressed in COS-7 cells with a variety of adrenergic ligands demonstrates a unique alpha 2AR pharmacology. Hybridization with somatic cell hybrids shows that the gene for this receptor is located on chromosome 2. Northern blot analysis of various rat tissues shows expression in liver and kidney. The unique pharmacology and tissue localization of this receptor suggest that this is an alpha 2AR subtype not previously identified by classical pharmacological or ligand binding approaches.

Cloning of the cDNA for the human beta 1-adrenergic receptor.

Screening of a human placenta lambda gt11 library has led to the isolation of the cDNA for the human beta 1-adrenergic receptor (beta 1AR). Used as the probe was the human genomic clone termed G-21. This clone, which contains an intronless gene for a putative receptor, was previously isolated by virtue of its cross hybridization with the human beta 2-adrenergic receptor (beta 2AR). The 2.4-kilobase cDNA for the human beta 1AR encodes a protein of 477 amino acid residues that is 69% homologous with the avian beta AR but only 54% homologous with the human beta 2AR. This suggests that the avian gene encoding beta AR and the human gene encoding beta 1AR evolved from a common ancestral gene. RNA blot analysis indicates a message of 2.5 kilobases in rat tissues, with a pattern of tissue distribution consistent with beta 1AR binding. This pattern is quite distinct from the pattern obtained when the beta 2AR cDNA is used as a probe. Expression of receptor protein in Xenopus laevis oocytes conveys adenylate cyclase responsiveness to catecholamines with a typical beta 1AR specificity. This contrasts with the typical beta 2 subtype specificity observed when the human beta 2AR cDNA is expressed in this system. Mammalian beta 1AR and beta 2AR are thus products of distinct genes, both of which are apparently related to the putative G-21 receptor.

Promoter rearrangements cause species-specific hepatic regulation of the glyoxylate reductase/hydroxypyruvate reductase gene by the peroxisome proliferator-activated receptor alpha.

In liver, the glyoxylate cycle contributes to two metabolic functions, urea and glucose synthesis. One of the key enzymes in this pathway is glyoxylate reductase/hydroxypyruvate reductase (GRHPR) whose dysfunction in human causes primary hyperoxaluria type 2, a disease resulting in oxalate accumulation and formation of kidney stones. In this study, we provide evidence for a transcriptional regulation by the peroxisome proliferator-activated receptor alpha (PPARalpha) of the mouse GRHPR gene in liver. Mice fed with a PPARalpha ligand or in which PPARalpha activity is enhanced by fasting increase their GRHPR gene expression via a peroxisome proliferator response element located in the promoter region of the gene. Consistent with these observations, mice deficient in PPARalpha present higher plasma levels of oxalate in comparison with their wild type counterparts. As expected, the administration of a PPARalpha ligand (Wy-14,643) reduces the plasma oxalate levels. Surprisingly, this effect is also observed in null mice, suggesting a PPARalpha-independent action of the compound. Despite a high degree of similarity between the transcribed region of the human and mouse GRHPR gene, the human promoter has been dramatically reorganized, which has resulted in a loss of PPARalpha regulation. Overall, these data indicate a species-specific regulation by PPARalpha of GRHPR, a key gene of the glyoxylate cycle.

The G0/G1 switch gene 2 is a novel PPAR target gene.

PPARs (peroxisome-proliferator-activated receptors) alpha, beta/delta and gamma are a group of transcription factors that are involved in numerous processes, including lipid metabolism and adipogenesis. By comparing liver mRNAs of wild-type and PPARalpha-null mice using microarrays, a novel putative target gene of PPARalpha, G0S2 (G0/G1 switch gene 2), was identified. Hepatic expression of G0S2 was up-regulated by fasting and by the PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. Surprisingly, the G0S2 mRNA level was highest in brown and white adipose tissue and was greatly up-regulated during mouse 3T3-L1 and human SGBS (Simpson-Golabi-Behmel syndrome) adipogenesis. Transactivation, gel shift and chromatin immunoprecipitation assays indicated that G0S2 is a direct PPARgamma and probable PPARalpha target gene with a functional PPRE (PPAR-responsive element) in its promoter. Up-regulation of G0S2 mRNA seemed to be specific for adipogenesis, and was not observed during osteogenesis or myogenesis. In 3T3-L1 fibroblasts, expression of G0S2 was associated with growth arrest, which is required for 3T3-L1 adipogenesis. Together, these data indicate that G0S2 is a novel target gene of PPARs that may be involved in adipocyte differentiation.

Structure of the glycosyl-phosphatidylinositol membrane anchor of the Leishmania major promastigote surface protease.

In common with many other plasma membrane glycoproteins of eukaryotic origin, the promastigote surface protease (PSP) of the protozoan parasite Leishmania contains a glycosyl-phosphatidylinositol (GPI) membrane anchor. The GPI anchor of Leishmania major PSP was purified following proteolysis of the PSP and analyzed by two-dimensional 1H-1H NMR, compositional and methylation linkage analyses, chemical and enzymatic modifications, and amino acid sequencing. From these results, the structure of the GPI-containing peptide was found to be Asp-Gly-Gly-Asn-ethanolamine-PO4-6Man alpha 1-6Man alpha 1-4GlcN alpha 1-6myo-inositol-1-PO4-(1-alkyl-2-acyl-glycerol). The glycan structure is identical to the conserved glycan core regions of the GPI anchor of Trypanosoma brucei variant surface glycoprotein and rat brain Thy-1 antigen, supporting the notion that this portion of GPIs are highly conserved. The phosphatidylinositol moiety of the PSP anchor is unusual, containing a fully saturated, unbranched 1-O-alkyl chain (mainly C24:0) and a mixture of fully saturated unbranched 2-O-acyl chains (C12:0, C14:0, C16:0, and C18:0). This lipid composition differs significantly from those of the GPIs of T. brucei variant surface glycoprotein and mammalian erythrocyte acetylcholinesterase but is similar to that of a family of glycosylated phosphoinositides found uniquely in Leishmania.

Cloning and characterization of the cDNA for the Brazilian Cratomorphus distinctus larval firefly luciferase: Similarities with European Lampyris noctiluca and Asiatic Pyrocoelia luciferases

Studies on firefly (Lampyridae) luciferases have focused on nearctic species of Photinus and Photuris and Euroasiatic species of Lampyris, Luciola, Hotaria, and Pyrocoelia. Despite accounting for the greatest diversity of fireflies in the world, no molecular studies have been carried out on the highly diverse genera from the neotropical region. Here we report the luciferase cDNA cloning for the larva of the Brazilian firefly Cratomorphus distinctus. The cDNA has 1978 bp and codes for a 547-residue-long polypeptide. Noteworthy, sequence comparison as well as functional properties show the highest degree of similarity with Lampyris noctiluca (93%) and Pyrocoelia spp. (91%) luciferases, suggesting a close phylogenetic relationship despite the geographical distance separating these species. The bioluminescence emission spectrum peaks at 550 nm and, as expected, is sensitive to pH, shifting to 605 nm at pH 6. The kinetic properties of the recombinant luciferase were similar to those of other firefly luciferases. © 2004 Elsevier Inc. All rights reserved.

Taxonomic and phylogenetic relationships between neotropical species of ticks from genus Amblyomma (Acari: Ixodidae) inferred from second internal transcribed spacer sequences of rDNA

The accurate specific identification of ticks is essential for the study, control and prevention of tick-borne diseases. Herein, we determined ribosomal nucleotide sequences of the second internal transcribed spacer (ITS2) of 15 Neotropical hard tick species of the genus Amblyomma Koch found in Brazil. Most of the studied ticks accidentally parasite humans and potentially act as vectors of zoonoses. Lengths of the ITS2 sequences ranged from 956 to 1,207 bp, whereas GC content varied from 62.4 to 66.9%. A matrix of ITS2 divergence was calculated with the ITS2 sequence data obtained showing divergence levels varying from 1.5 to 28.8%. The analysis indicated that this molecular marker can be useful for Amblyomma-specific identification. Phylogenetic inferences based on the ITS2 sequences were used to assess some issues in subgenus taxonomy. © 2007 Entomological Society of America.